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1.0  ABSTRACT 


A method  has  been  developed  to  provide  improved  trajectory  calculations 
for  helicopter- launched  missiles  based  on  an  earlier  U.S.  Amy  trajectory 
program.  In  the  earlier  method  the  helicopter  rotor  flow  field  is  simulated 
but,  at  a given  instant  of  time,  the  flow  is  assumed  to  be  uniform  over  the 
whole  missile.  By  incorporating  a three-dimensional  panel  method  the  forces 
and  moments  due  to  the  flow-field  nonuniformity  are  calculated  and  applied  in 
the  form  of  a perturbation  about  the  uniform  flow  values. 

The  uniform  flow  values  can  be  estimated  by  using  experimental  results 
incorporated  into  the  basic  trajectory  calculation  and  so  in  this  way  the  pre- 
dominant effects  of  the  flow-field  separation  can  be  accurately  accounted  for. 

A relatively  crude  model  of  the  flow  separation  for  the  perturbation  solution 
is  therefore  used  in  which  a blowing  velocity  is  imposed  over  the  part  of  the 
missile  for  which  the  flow  would  be  separated. 

This  method  has  been  applied  to  the  calculation  of  the  trajectory  of  the 
2.75-Inch  Rocket  launched  from  a hovering  helicopter  using  two  different  down- 
wash  distributions.  It  is  found  that  the  modification  to  the  trajectory  intro- 
duced by  the  calculated  effects  due  to  the  flow-field  nonuniformity  is  relatively 
small  but  it  can  still  be  significant  as  regards  accuracy  at  the  target. 
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3.0  INTRODUCTION 


The  basic  aim  of  the  work  described  in  this  report  is  to  obtain  an 
improved  trajectory  calculation  for  a finned  missile  launched  from  a heli- 
copter in  flight.  In  particular,  the  study  is  aimed  at  determining  the  effect 
on  the  trajectory  of  the  nonuniform  character  of  the  velocity  field  of  the 
rotor  that  is  experienced  by  the  missile  at  any  instant  of  time.  Existing 
procedures  (references  1 and  2)  account  for  the  fact  that  the  flow  environ- 
ment changes  as  the  missile  advances  radially  outward  through  the  rotor  field. 
However,  these  procedures  assume  either  that  the  velocity  field  experienced  by 
the  missile  at  any  particular  instant  of  time  is  uniform  or  that  a simple 
correction  to  a uniform  flow  can  be  made  (reference  2).  The  basic  calcula- 
tional  tool  used  to  evaluate  the  effects  of  flow  nonuniformity  is  the  panel 
method  of  reference  3,  which  is  described  briefly  in  Section  5.0. 

It  is  clear  that  during  the  launch  phase  a missile  in  the  strong  down- 
wash  field  of  a rotor  experiences  large  angles  of  attack  and  that  the  flow 
about  the  missile  is  largely  separated.  Accordingly,  the  early  stages  of  the 
present  work  concentrated  on  developing  a potential -flow  model  of  separated 
flow  that  could  be  implemented  by  means  of  a panel  method.  Unfortunately,  the 
attempt  was  unsuccessful.  This  was  partly  due  to  the  very  complicated  nature 
of  three-dimensional  separation  (reference  4)  and  partly  due  to  uncertainties 
of  constructing  potential-flow  models  of  separation  (reference  5).  In  any 
event  the  task  is  clearly  beyond  the  scope  of  the  present  work. 

An  additional  consideration  is  the  rather  good  accuracy  of  the  results 
obtained  by  existing  methods.  In  the  near- to-medi urn  future  no  method  basf.d 
on  fundamental  flow  analysis  will  be  able  to  match  the  accuracies  obtained 
by  the  more  practically-oriented  techniques  of  references  1 and  2.  Moreover, 
the  fundamental  procedure  would  require  orders  of  magnitude  more  computing 
time.  Once  it  had  been  decided  that  the  present  work  would  be  directed 
towards  development  of  a practical  method,  i.e.,  one  that  offers  a reasonable 
hope  of  improving  on  existing  methods  with  an  acceptable  computing  time,  the 
permissible  outline  of  the  method  became  fairly  clear. 
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4.0  GENERAL  DESCRIPTION  OF  THE  METHOD 
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L. 


The  essential  property  of  the  method  of  this  report  is  that  it  calculates 
a perturbation  to  a uniform  onset-flow  solution,  which  thus  becomes  the  total 
solution  as  the  flow  nonuniformity  vanishes.  Specifically,  it  is  assumed  that 
tabular  data  exists  that  permit  determination  of  force  and  moment  coefficients 
due  to  a uniform  freestream  at  any  inclination  to  the  missile  axis.  The 
source  of  this  data  is  not  considered  here.  Ideally,  it  would  be  gathered 
from  wind  tunnel  tests.  The  method  also  relies  on  existing  programs  of  the 
type  of  references  1 and  2,  which  not  only  provide  machinery  for  integrating 
the  missile's  trajectory  given  its  force  and  moment  history,  but  also  contain 
subroutines  that  give  the  downwash  field  of  the  rotor  as  a function  of  radial 
location  and  relative  blade  position.  Essentially,  the  panel  method  is  incor- 
porated into  this  last  type  of  program  where  it  computes  perturbations  to  the 
uni  form- flow  forces  and  moments. 

At  any  given  time  the  known  missile  position  is  used  together  with  the 
rotor-field  subroutine  to  give  the  onset  flow  at  all  points  of  the  missile. 

One  point  is  selected  where  the  magnitude  and  direction  of  the  onset  flow 
velocity  is  felt  to  be  particularly  important.  This  velocity  is  used  with 
the  tabular  data  to  obtain  uni  form- flow  forces  and  moments.  The  panel  method 
calculates  perturbations  to  the  forces  and  moments  based  on  differences 
between  the  nonuniform  onset  flow  due  to  the  rotor  and  the  above  uniform 
onset  flow,  as  described  in  Section  6.0.  A somewhat  open  question  is  the 
selection  of  the  point  where  the  uniform  onset-flow  velocity  is  evaluated. 

For  a missile  with  rear-mounted  fins,  it  appears  that  the  principal  effect  of 
the  downwash  field  on  the  flight  path  is  the  pitch  up  resulting  from  the  down 
force  on  the  fins.  Accordingly,  the  fin  location  has  been  selected  as  this 
velocity  reference  point.  Some  calculated  results  showing  the  importance 
of  this  choice  are  presented  in  Section  7.0. 

The  above  formulation  has  several  important  advantages.  The  most 
important  is  that  the  principal  effects  of  the  calculationally  very  difficult 
phenomenum  of  separation  are  contained  in  the  uniform-flow  solution,  which  can 
be  obtained  routinely  in  a wind  tunnel.  Thus  the  calculational  effort  is  con- 
fined to  the  nonuniform  field  effect,  which  is  beyond  the  capability  of  a 
wind  tunnel.  Since  it  is  only  a perturbation,  the  calculated  flow  can  use  a 
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very  simple  separation  model,  similar  to  that  of  reference  6,  without 
introducing  significant  inaccuracies.  Finally,  the  panel-method  calculation 
need  be  carried  out  only  while  the  missile  is  in  the  field  of  the  rotor,  which 
represents  a very  small  fraction  of  the  total  flight  time.  Once  the  missile 
has  left  the  rotor  field  the  trajectory  calculation  becomes  very  simple.  A 
more  fundamental  attempt  at  flow  modelling  would  require  use  of  the  panel 
method  over  the  entire  flight  path  with  a consequent  massive  increase  in 
computer  time. 
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5.0  THE  PANEL  METHOD 


The  evaluation  of  the  nonunifonn-flow-field  effects  involves  only  the 
calculation  of  small  perturbations  to  the  uniform-flow  forces  and  moments 
induced  by  the  nonuniform  part  of  the  flow  field.  As  stated  in  the  previous 
section,  the  velocity  reference  point  at  which  the  uni  form- flow  velocity  is 
evaluated  is  taken  to  be  at  the  missile  fin  location.  Therefore,  the  non- 
uniform  part  of  the  flow  field  will  vanish  at  this  point  and  so,  to  a good 
approximation,  the  flow  about  the  missile  due  to  the  nonuniform  component  of 
the  flow  field  will  be  unaffected  by  the  presence  of  the  fins.  Accordingly, 
it  was  decided  to  apply  the  panel-method  calculation  to  the  missile  without 
its  fins.  This,  in  turn,  permits  the  use  of  the  three-dimensional  nonlifting 
panel  method  described  in  reference  3,  since  no  lifting  surfaces  remain  on 
the  missile.  This  method  is  simpler,  more  flexible  and  faster  than  a lifting 
method  whose  use  would  add  significantly  to  the  complexity  and  computing  time 
of  the  resulting  computer  program. 

For  full  details  of  the  panel  method  used,  the  reader  is  referred  to 
reference  7.  Only  a brief  description  will  be  given  here,  together  with 
details  of  the  modifications  to  the  existing  computer  program  which  were 
necessary  to  produce  a method  suitable  for  the  current  application.  The  model 
used  represents  the  missile  by  a number  of  flat  quadrilateral  panels,  mostly 
nearly  rectangular  (see  fig.  1),  on  each  of  which  there  is  a constant  source 
strength.  A matrix  of  influence  coefficients,  representing  the  normal  velocity 
induced  on  all  of  these  panels  by  the  singularity  distribution  on  each  of  the 
panels,  is  then  calculated.  This  matrix  is  then  used  to  determine  the  source 
strengths  in  such  a way  that  the  boundary  condition  of  prescribed  total  normal 
velocity  is  satisfied  on  each  of  the  panels.  Additional  matrices  representing 
the  tangential  velocities  induced  by  each  panel  are  calculated  so  that,  given 
the  singularity  distribution,  the  entire  velocity  distribution  over  the  body 
can  be  computed.  This  calculation  can  be  carried  out  routinely  regardless  of 
the  nature  of  the  flow  field  incident  to  the  body. 

In  order  to  reduce  the  number  of  panels  which  need  to  be  considered, 
advantage  is  taken  of  the  symmetry  about  the  vertical  lengthwise  plane  which 
the  missile  possesses.  (Clearly  the  missile  has  a good  deal  more  symmetry 


degree.)  The  flow  field  itself  is  not  symmetric  about  this  plane  but  it 
can  be  separated  into  two  components,  one  symmetric  and  the  other  antisym- 
metric whose  sum  is  equal  to  the  physical  flow  field.  Separate  influence 
coefficients  are  required  for  these  two  components,  but  their  solution  requires 
less  computer  time  than  would  be  required  to  solve  one  matrix  equation  with 
twice  as  many  unknowns, which  would  be  required  without  the  use  of  this  sym- 
metry option. 

The  calculation  of  the  missile  behavior  through  the  helicopter  rotor 
field  requires  a panel  method  calculation  at  each  successive  point  of  the 
trajectory.  Since  the  normal  and  tangential  velocity  matrices  are  dependent 
only  on  the  body  geometry,  it  is  not  necessary  to  recompute  them  for  each  time 
step.  In  addition,  a considerable  amount  of  computing  time  can  be  saved  by 
storing  the  inverse  normal  velocity  matrices  for  use  at  each  successive  time 
step.  Thus,  the  amount  of  calculation  required  at  each  time  step  is  a small 
fraction  of  that  required  for  the  initial  panel  method  calculation. 

The  basic  panel  method  was,  therefore,  modified  to  permit  the  tangential 
velocity  matrices  and  the  inverse  normal  velocity  matrices  to  be  stored  on 
disk.  By  greatly  restricting  the  maximum  number  of  panels  which  could  be 
used,  as  compared  with  the  basic  method,  (see  Section  7.0)  it  was  possible  to 
enable  one  entire  influence  matrix  to  be  held  in  the  high-speed  core  at  any 
given  time.  This  permits  the  use  of  an  in-core  matrix  solution  routine  which 
leads  to  a further  reduction  in  the  computing  time  by  reducing  the  number  of 
disk  transfers  required. 

In  addition  to  these  modifications,  the  program  was  further  simplified 
to  remove  a number  of  the  external  storage  units.  The  basic  program  used 
these  units  for  transferring  data  from  one  segment  of  the  program  to  another, 
but  with  the  reduced  panel  number  used  here,  these  data  can  now  be  retained 
in  the  core. 


6.0  CALCULATION  OF  THE  PERTURBATION  FORCES 


As  stated  in  the  Introduction,  the  purpose  of  the  panel -method  calcula- 
tion is  to  provide  perturbations  to  the  aerodynamic  forces  and  moments  due  to 
the  flow  nonuniformities  acting  on  the  missile.  Therefore,  it  is  necessary 

to  calculate  the  pressure  distribution  due  to  both  the  total  onset  flow,  Uj, 

-►  * 
and  the  uniform  onset  flow,  u^,  which  is  defined  by  evaluating  the  velocity 

at  a predetermined  reference  point  on  the  body.  Typically,  the  missile  is  in 

the  rotor  downwash  for  about  0.1  seconds  and  the  principal  effect  of  this 

velocity  field  on  a missile  with  rear  mounted  fins  is  to  induce  a pitching 

moment  caused  by  the  aerodynamic  forces  on  the  fins.  This  velocity  reference 

point  is,  therefore,  taken  to  be  near  the  tail. 

A simplified  model  of  the  separated -flow  effect  is  used  in  which  a 
positive  normal  blowing  velocity  is  introduced  over  panels  for  which  the  flow 
is  separated.  The  separated-flow  region  is  simply  defined  by  that  region  for 
which  the  normal  component  of  the  onset  velocity  is  positive,  and  the  magni- 
tude of  the  blowing  velocity  is  taken  to  be  equal  to  the  normal  component  of 
the  onset  velocity.  When  defined  in  this  way  the  blowing  velocity  tends  to 
zero  at  the  edges  of  the  separated-flow  region  and,  if  applied  to  two- 
dimensional  flow  past  a cylinder,  it  would  result  in  a wake  behind  the  cylinder 
whose  downstream  thickness  tends  to  the  diameter  of  the  cylinder.  The  precise 
definition  of  these  separated  flow  parameters  can  be  modified,  but  the  results 
presented  in  the  next  section  indicate  that  there  is  not  a great  deal  of 
dependence  on  the  separation  modelling. 

Two  blowing  velocities  are  therefore  defined  for  each  panel  of  the  body, 
by  and  b^,  corresponding  to  the  total  and  the  uniform  onset  flows, 
respectively.  In  the  attached-flow  regions  these  blowing  velocities  are 
defined  to  be  zero.  To  achieve  a solution  with  these  blowing  distributions,  the 
right-hand  sides  for  the  normal  velocity  equations  are  defined  by 

-Uj  • n + bj  for  the  total  flow 

and  (6.1) 

-u^  • n + b^  for  the  uniform  flow 

0 0 

A 

where  n is  the  unit  normal  vector. 
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The  two  vectors  defined  by  evaluating  (6.1)  for  each  panel  on  the 
missile  are  each  separated  into  two  components,  a symmetrical  and  an  anti- 
symmetrical  component.  The  respective  source  densities  required  to  generate 
these  normal  velocities  over  the  body  are  then  calculated  using  the  inverse 
normal  velocity  matrices  and  the  corresponding  velocity  distributions  are 
computed  using  the  tangential  velocity  matrices.  At  this  stage  the  symmetric 
and  antisymmetric  components  are  recombined  to  give  the  surface  velocity 
distributions  Vj  and  v^  over  the  entire  missile. 

By  making  the  assumption  that  the  freestream  dynamic  pressure  for  any 
streamline  impinging  on  the  missile  at  any  instant  of  time  is  constant  and 
equal  to  the  uniform-flow  value,  then  the  pressure  coefficient  in  the 
attached  flow  region  can  be  written  as 


0 


where  v^  is  the  magnitude  of  the  velocity  at  the  velocity  reference  point. 
This  formula  cannot,  however,  be  used  in  the  separated -flow  region  since  an 
artificial  blowing  velocity  has  been  introduced  to  model  the  flow  outside 
the  separated -flow  region.  In  this  region,  at  a given  lengthwise  station, 
the  pressure  is  taken  to  be  constant  with  its  value  given  by  the  average  of 
the  values  at  the  control  points  on  either  side  of  the  separated  region.  The 
computer  program  is,  therefore,  set  up  to  calculate  the  pressure  coefficients 
for  both  the  uniform  and  the  total  flows  based  on  equation  (6.2)  after  which 
the  separated-flow  values  are  overwritten  as  outlined  above. 

The  perturbation  forces  and  moments  are  then  calculated  by  integrating 
the  difference  between  the  total-flow  pressure  and  the  uniform-flow  pressure 
over  the  whole  missile.  These  are  then  added  to  the  corresponding  uniform- 
flow  forces  and  moments  (obtained  experimentally  or  otherwise)  before  pro- 
ceeding with  the  trajectory  calculation. 


7.0  DISCUSSION  OF  CALCULATED  RESULTS 


A measure  of  the  normal -plane  impact  error  in  a missile  trajectory  is 
provided  by  the  errors  in  the  flight-path  angle  at  the  end  of  the  boosted 
stage  of  the  flight.  Thus  the  results  presented  here  cover  only  this  seg- 
ment of  the  missile  flight  with  both  instantaneous  pitch  angle  and  the  flight- 
path  angle  in  the  vertical  plane  being  plotted.  The  principal  effect  of  the 
imposed  downwash  field  occurs  in  the  vertical  plane  and  so  the  horizontal 
deviation  has  not  been  plotted  here  although  these  data  are  calculated  by  the 
program. 

ilesults  for  two  downwash  fields  will  be  presented  in  detail:  the  first 
a linear  downwash  in  which  the  vertical  velocity  is  proportional  to  the  dist- 
ance from  the  axis  of  the  rotor  and  the  second  a simulated  helicopter  rotor 
field  based  on  a helical  vortex  model. 

7.1  Linear  Downwash 

The  linear  downwash,  which  is  zero  on  the  rotor  axis  and  116  ft/sec  at 
a distance  of  20  feet  from  the  axis,  was  considered  by  Figueras  and  Bauman 
(ref.  1)  in  an  earlier  study  of  the  effect  of  a downwash  distribution  on  a 
helicopter-launched  missile.  This  case,  therefore,  provides  a useful  compar- 
ison with  the  earlier  results. 


The  missile-launch  point  is  taken  to  be  on  the  helicopter  centerline  at 
a distance  of  3 feet  ahead  of  the  rotor  axis,  and  figures  2 and  3 show  the 
instantaneous  missile  pitch  and  flight-path  angle  respectively,  where  three 
curves  have  been  plotted  in  each  figure.  The  first  is  that  obtained  from  the 
basic  trajectory  program  in  which  the  downwash  field  is  evaluated  at  the 
missile  center  of  gravity  and  is  assumed  to  apply  uniformly  over  the  whole 
length  of  the  missile.  The  second  curve  represents  a minor  modification  to 
this  program  in  which  the  velocity  is  still  assumed  to  be  uniform  over  the 
missile,  but  the  reference  point  at  which  it  is  evaluated  is  taken  to  be  near 
the  missile  tail.  The  third  curve  uses  this  same  reference  velocity  but  it 
also  includes  the  panel -method  evaluation  of  the  effects  of  the  flow  field 
nonuniformity.  It  can  be  seen  from  these  figures  that  a significant  correc- 
tion in  the  missile  trajectory  is  achieved  simply  by  moving  the  velocity  ref- 
erence point  from  the  missile  center  of  gravity  to  the  tail.  A further 
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correction  is  then  achieved  by  allowing  for  the  flow-field  nonuni formlty. 

For  this  particular  downwash  distribution,  both  of  these  corrections  act  In 
the  same  direction  producing  a nose-down  pitch  relative  to  the  basic  method. 
This  result  Is  clearly  what  would  be  expected  for  this  simple  downwash  field 
since  the  evaluation  of  the  uniform  velocity  at  the  tall  gives  a reduced 
vertical  component  and  hence  a reduced  pitch  up  due  to  the  Influence  of  this 
velocity  on  the  missile  fin,  and  the  nonunlform  flow  effect  contributes  further 
nose-down  pitching  moment.  A total  correction  of  about  1.3  degrees  (22  mils) 

In  the  flight-path  angle  at  the  end  of  the  boosted  flight  Is  obtained  between 
the  basic  method  and  the  panel-method  calculation. 

It  should  be  noted  that  there  Is  a difference  between  the  results  pre- 
sented here  for  the  basic  method  and  those  presented  in  fig.  10  of  reference  1. 
The  results  have  the  same  qualitative  behavior  but  both  the  magnitude  and 
frequency  of  the  oscillations  occurring  in  the  pitch  angle  (fig.  2)  differ. 
Despite,  as  far  as  possible,  using  the  same  missile  characteristics  and 
launch  point.  It  was  not  possible  to  reproduce  exactly  these  earlier  results. 
Therefore,  It  was  concluded  that  this  discrepancy  must  be  due  to  the  use  of 
slightly  different  dynamic  or  aerodynamic  data  within  the  trajectory 
calculation. 

For  this  particular  downwash  the  calculation  appears  to  be  Insensitive 
to  the  separated -flow  model  which  Is  used  since  the  results  are  virtually 
unchanged  when  the  flow  Is  assumed  to  be  totally  attached.  In  addition,  the 
calculation  is  relatively  Insensitive  to  the  number  of  panels  used  to  repre- 
sent the  missile.  The  results  presented  here  were  all  obtained  using  60 
panels  over  half  of  the  missile,  10  In  the  axial  direction  and  6 In  the  cir- 
cumferential direction.  By  doubling  the  number  of  circumferential  panels, 
but  using  the  same  axial-panel  distribution,  the  results  for  the  triangular 
downwash  distribution  of  figs.  2 and  3 were  graphically  Identical. 

L 

7.2  Helicopter  Rotor  Field 

A more  realistic  model  of  the  downwash  below  a helicopter  rotor  Is 
^ provided  by  a set  of  routines  used  In  the  U.S.  Army  6-degree-of- freedom 

I trajectory  program  used  to  simulate  the  2. 75-Inch  rocket.  This  model  uses 

I a helical  vortex  which  Is  assumed  to  leave  each  rotor  blade.  The  location 
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and  strength  of  these  vortices  is  given  by  a set  of  data,  while  the  program 
calculates  the  resulting  velocity  field  including  modifications  due  to  the 
presence  of  the  helicopter  fuselage.  These  routines  are  used  to  evaluate  the 
flow  field  at  a number  of  points  on  the  missile  axis,  these  values  are  then 
used  to  compute  the  normal  velocity  distribution  on  the  missile  surface. 

Figures  4 and  5 show  the  pitch  and  flight-path  angle  for  the  2.75-inch 
missile  launched  through  this  downwash  field  with  no  initial  incidence.  Again 
these  figures  show  the  curves  obtained  for  the  uniform-flow  calculations  with 
the  reference  velocity  evaluated  both  at  the  missile  center  of  gravity  and  at 
the  missile  tail  in  addition  to  the  nonuniform  panel -method  solution.  For 
this  downwash  distribution  it  can  be  seen  that  moving  the  velocity  reference 
point  alone  overcompensates  for  the  error  in  the  basic  method  as  compared  with 
the  panel-method  solution.  The  correction  introduced  by  the  flow-field  non- 
uniformity for  this  case  amounts  to  about  0.5  degrees  or  9 mils  at  the  end  of 
the  boosted  flight. 

As  with  the  triangular  velocity  distribution,  there  is  no  noti cable 
difference  in  these  figures  if  the  panel  number  is  increased  to  120  and  so  it 
was  concluded  that  60  panels  provide  an  adequate  missile  representation.  In 
this  case,  however,  there  is  some  dependence  on  the  separation  algorithm  as 
can  be  seen  from  figures  6 and  7 which  show  the  results  for  which  the  flow  is 
assumed  to  be  attached.  In  this  attached  flow  case  there  is  virtually  no 
difference  between  the  results  with  and  without  the  panel  method  provided 
that  the  tail  is  used  as  the  velocity  reference  point. 

Since  the  helicopter  rotor  field  is  dependent  on  the  position  of  the 
rotor  blades,  the  panel  method  has  been  applied  with  various  rotor  positions. 
Figures  8-11  show  the  results  obtained  for  rotor  positions  varying  by  45° 
and  90°  from  the  position  used  for  the  previous  results.  For  these  particular 
combinations  of  launch  conditions  and  rotor  positions,  there  is  only  a small 
change  in  the  calculated  missile  trajectory.  However,  due  to  the  highly 
nonuniform  character  of  the  flow  field,  there  could  be  cases  for  which  the 
rotor  position  can  have  a significant  influence.  Similarly,  the  initial 
launch-angle  can  have  a significant  effect  on  which  part  of  the  flow  field  the 
missile  encounters.  Two  figures,  12  and  13,  plot  the  pitch  and  flight-path 
angles  for  a missile  launched  at  an  incidence  of  7°.  In  this  case  it  can  be 
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seen  that  a nose-down  pitch  correction  Is  provided  by  the  panel  method 
whereas  for  the  horizontally  launched  missile  the  panel  method  provides  a 
nose-up  pitch  correction. 


All  of  the  results  presented  so  far  have  been  derived  using  a time  step 
of  0.025  in  the  trajectory  calculation  which,  while  providing  a useful  qual- 
itative comparison  between  the  different  cases  and  methods.  Is  too  large  to 
permit  reliable  results  to  be  obtained.  Some  results  for  smaller  time  steps 
are  presented  In  figures  14-17.  It  can  be  seen  that,  by  reducing  the  time 
step  to  0.0125  seconds,  both  the  Initial  peak  In  the  nose-up  pitch,  which  occurs 
at  about  0.27  seconds,  and  the  subsequent  oscillations  In  the  missile  attitude 
are  Increased  In  magnitude.  Similarly,  the  flight-path  angle  (figure  15)  Is 
also  Increased  over  the  whole  flight  time  plotted  here  for  each  of  the  three 
methods  with  the  largest  Increase  occurring  In  the  basic  solution.  Figures 
16  and  17  show  that  a further  reduction  In  the  step  size  to  0.005  produces  a 
less  significant  change  In  the  solution.  It  can,  therefore,  be  concluded 
that  for  most  cases  a time  step  of  about  0.0125  seconds  could  be  used.  For 
greater  accuracy  the  time  step  could  be  further  reduced  to  about  0.005 
seconds  although  this  will  lead  to  a corresponding  Increase  In  computer  time. 

For  example,  the  total  computing  time  required  on  an  IBM  370/3033  system  for 
a time  step  of  0.005  seconds  was  about  100  seconds  of  which  about  60  seconds 
was  required  for  all  of  the  helicopter  rotor  field  evaluations  and  about  30 
seconds  was  required  for  all  of  the  panel -method  calculations.  The  correspond- 
ing time  for  the  basic  trajectory  program  with  the  same  time  step  was  about 
10  seconds. 


8.0  PROGRAM  INPUT  INSTRUCTIONS 

The  input  Instructions  for  the  modified  version  of  the  trajectory  program 
can  be  conveniently  divided  into  two  parts,  the  first  providing  control,  flags 
and  input  data  for  the  panel  method  and  the  second  providing  the  data  required 
for  the  trajectory  calculation.  The  latter  is  very  similar  to  the  data  required 
for  the  basic  trajectory  program  with  only  a few  changes  and  additions  being 
necessary,  but  a brief  description  of  both  data  sets  will  be  provided  here. 

8.1  Panel  Method  Input 

The  panel  method  is  executed  once  at  the  start  of  the  program  in  order 
to  read  in  the  missile  geometry  data  and  calculate  the  normal  and  tangential 
velocity  matrices.  At  this  time  the  normal  velocity  matrices  are  factorized 
into  upper  and  lower  triangular  matrices  to  permit  the  rapid  calculation  of 
the  solutions  for  additional  right-hand  sides  required  for  each  time  step. 

The  tangential  velocity  matrices  and  the  factorized  normal  velocity  matrices 
are  then  stored  on  external  units  for  use  as  the  calculation  proceeds.  Since 
these  quantities  are  dependent  solely  on  the  missile  panelling  which  will,  in 
general,  remain  unchanged  from  one  run  to  the  next,  an  option  is  provided 
for  bypassing  much  of  this  initial  calculation  by  using  the  matrices  stored 
during  a previous  run. 

As  mentioned  in  Section  5.0,  the  panel  method  makes  the  basic  assumption 
that  the  missile  is  symmetric  about  a plane  through  the  missile  axis.  In 
fact,  as  far  as  the  geometric  input  is  concerned,  the  missile  is  assumed  to 
be  axisymmetric  so  that,  if  the  missile  profile  is  specified,  the  program 
generates  the  circumferential  panel  distribution.  The  array  dimensions 
currently  defined  allow  for  up  to  10  panels  in  the  axial  direction  with  up 
to  18  panels  in  the  circumferential  direction,  over  half  of  the  missile. 

Thus  up  to  eleven  pairs  of  coordinates  must  be  input  to  define  the  axial  panel 
distribution. 

The  data  for  the  panel  method  is  input  on  unit  5 and  the  cards  required 


Card  1 — Control  Card 


Format  (212) 

cc 

Variables 

Remarks 

1-2 

MODEl 

=0 

for  an  initial  run  with  new  panel 
distribution 

=1 

for  subsequent  runs  using  geometry  data 
stored  on  units  8,  9 and  13  during  a 

previous  run 

3-4 

M0DE2 

=0 

if  the  helicopter  downwash  is  to  be  used 

=1 

if  a linear  downwash  is  to  be  used 

Card  2 - Downwash  Parameters 
This  card  in  input  only  if  MODEZ^O 
Format  (2F10.6) 

cc  Variables  Remarks 

1-10  DNl  A linear  downwash  velocity  is  defined  by 

11-20  DN2  DN1+DN2*DDISS  where  DDISS  is  the  hori- 

zontal distance  from  the  rotor  axis 

The  remaining  cards  must  be  input  only  if  M0DE1=0. 

Card  3 - Title  Card 

Format  C15A4,A4) 

cc  Var^  V "es  Remarks 

1-60  HEDR  Alpha-numeric  title  and  case  number  for 

61-64  KASE  printing  at  head  of  panel  method  output 

Card  4 - Input  Control 

Format  (2I5,2F10.5) 


cc 

Variables 

Remarks 

1-5 

N 

Number  of  axial  panels  (N  £ 10) 

6-10 

M 

Number  of  circumferential  panels  over  half 
of  the  missile  (M  £ 18) 

11-20 

B=1.0 

Length  of  axes  defining  elliptic  missile 

21-30 

C=-1.0 

cross-section  (N.B.  C must  be  negative 
to  ensure  that  the  missile  is  panelled  in 
the  correct  order.) 
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Cards  5-7  - X-Z  Input 
Format  (6F10.0) 


r 

I 


The  input  consists  of  coordinate  pairs  (X,,  Z, , X.,  Zp,  Xj.  , , Z,.  ,) 
defining  the  profile  of  the  body  to  be  generated  by  the  prografflT  ThKi 
points  define  the  axial  panel  distribution. 

8.2  Trajectory  Input 

The  data  required  for  the  trajectory  calculation  is  input  on  unit  3 in 
Namelist  Format.  A sample  input  list  will  be  given  together  with  a brief 
explanation  of  the  variables: 

&PUTT 


ALTIN  = -100.0, 

Initial  launch  altitude 

CKO  = 0.0, 

Thrust  misalignment 

CKl  = 0.0, 

Vectors 

CK4  = 6360.0, 

Related  to  missile  thrust 

DELgi2  = 0.00625, 

Half  time  step 

DELTT  =0.0125, 

Time  step 

DIVANG  =0.0, 

Dive  angle  of  helicopter  (degrees) 

DNSCH  = 0.5, 

Constant  in  dynamic  pressure  equation 

DPIT  = 0.0, 

DYAW  = 0.0, 

Thrust  misalignment  angles  (radians) 

DXLN  = 3.0, 

DYLN  = 0.0, 

DZLN  =7.0, 

Missile  launch  position  relative  to  center  of 
helicopter  rotor 

DZZ  = 0.0, 

D2RFF  = 0.2292, 

Missile  diameter  in  feet 

FACT0R  =0.0, 

FAC  = -6.56, 

Used  to  calculate  trim  pitch  attitude  of  helicopter 

FG  = 0.0, 

FT0RM  = 0.0, 

Specifies  units  (0.0  for  English  units  or  1.0  for 
metric  units) 

GAIN!  = 0.0, 


Not  used 


6AIN2  =0.0,  Not  used 

GREFF  = 32.173996,  Gravitational  constant  at  sea  level 

HELIC  Control  flag: 

= 0.0,  for  no  downwash  (and  no  panel  method  calculation) 

= 1.0,  for  calculation  of  downwash  effects  using  panel 

methods 

= -1.0,  for  calculation  of  downwash  effects  using  basic 

trajectory  program 

= 2.0,  for  calculation  of  downwash  effects  using  basic 

trajectory  program  but  with  reference  velocity 

evaluated  at  the  missile  tail 


I INPR0  = 0, 

i 

I 

] 

I INPUT0  = 1 , 


INT0PT  = 1 , 
ITUMB  = 1, 


JET  = 0, 
N = 25, 


NN  = 1, 

PHIZD  = 0.0, 
PSIZD  = 0.0, 
THTZD  = 0.0, 

PITMAL  =0.0, 
YAWMAL  = 0.0, 

QELNCH  =0.0, 

RAMP  = 0.0, 

WCF  = 0.0, 

WDF  = 0.0, 

REF  = 20902704.0 

R0TRAD  =22.0 

SCARF  = 0.00219 

SPNRAT  =0.0 


Indicates  that  missile  reference  axes  correspond 
to  missile  principal  axes 

Restart  parameter,  set  = 1 for  new  case 

Not  used 

Specifies  axes  fixed  with  respect  to  rotating 
missile 

Specifies  no  jet  damping 

Number  of  variables  to  be  integrated  by  Runge-Kutta 
integration  (leave  fixed) 

First  variable  to  be  integrated  (leave  fixed) 

Initial  missile  launch  angles  (degrees) 


Pitch  aid  yaw  tip  rates  as  missile  leaves  launcher 
(radians/sec) 

Missile  launcher  angle  (degrees) 

Horizontal  wind  components  defined  by 
WCF+RAMP*TIME  in  the  cross-range  direction  and 
WDF  down  range 


Radius  of  earth  (feet)' 


SUS0N  =0.0 
SUS0FF  = 1.551 
TIME  = 0.092, 
TFINAL  = 1.592, 
TR0T  =0.092, 

TLL  = 0.092, 
TPRINT  = 0.025, 
TPL0T  =0.025, 

T1  = 0.22, 

VINIT  = 0.0, 

XIND0T  = 0.0, 
YIND0T  = 0.0, 
ZIND0T  = 0.0, 

XLENG  = 4.4, 

XMASSO  = 0.61727, 

XMAX  = 0.0, 

XMIN  = 0.0, 

XMV  = 109.0, 

XTAR  = 13000.0, 
YTAR  = 0.0, 

ZTAR  = 0.0, 

YJET  = 0.0, 

ZJET  = 0.0, 

ZLL  = 0.0, 

ZLM  = 0.0, 

ZLN  = 0.0, 

&END 


Booster  motor  burn 

Times  (seconds) 

Missile  launch  time  (seconds) 

End  of  calculation  (seconds) 

Used  to  vary  rotor  blade  phase  at  time  of  missile 
launch 

Time  missile  leaves  launcher  (seconds) 

Time  interval  at  which  trajectory  data  is  printed 

Time  interval  at  which  missile  pitch  and  flight  path 
angles  are  stored  for  subsequent  plotting 

Time  at  which  panel  method  calculations  are  discon- 
tinued (seconds) 

Helicopter  velocity  at  missile  launch  time  (ft/sec) 


Launcher  length  (feet) 

Initial  missile  mass 
Not  used 
Not  used 

Missile  velocity  out  of  launch  tube  (ft/sec) 

Target  range, 
cross-range  and 
altitude  (feet) 

Jet  damping  terms 
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On  completion  of  one  case  the  program  begins  another  calculation  with 
the  same  missile  geometry.  Execution  is  terminated  by  specifying  a negative 
value  for  TIME,  and  so  the  following  cards  must  be  included  on  the  end  of 
this  data  set. 

&pun 

TIME  = -1.0, 

&END 

8.3  External  Units 

This  program  requires  7 external  units  whose  unit  numbers  and  contents 
are  listed  below: 

1.  Unit  2:  Used  for  storing  trajectory  data  for  subsequent  analysis 
or  plotting. 

2.  Unit  4:  Temporary  storage  unit  used  within  the  panel  method. 

3.  Unit  8:  Used  to  store  the  inverse  normal  velocity  matrices. 

4.  Unit  9:  Used  to  store  geometric  data  calculated  during  the  first 
call  to  the  panel  method. 

5.  Unit  10:  Helicopter  rotor  downwash  data  to  be  input  to  the 
trajectory  program. 

6.  Unit  11:  Used  to  store  pitch  and  flight- path  angle  data  for  sub- 
sequent analysis  or  plotting. 

7.  Unit  13:  Used  to  store  tangential  velocity  matrices. 

As  mentioned  in  subsection  8.1,  the  output  from  units  8,  9 and  13  can  be 
used  for  subsequent  runs  thereby  avoiding  the  need  to  recompute  the  influence 
matrices  and  their  inverses  if  they  remain  unchanged  from  one  run  to  the  next. 

8.4  Changes  Necessary  to  Run  an  Alternative  Missile 

The  panel-method  calculation  is  set  up  to  handle  any  axisymmetric  missile 
geometry,  but  the  trajectory  program  with  which  it  is  linked  is  based  on  the  U.S. 
Am\y  2. 75- inch  rocket.  Therefore,  for  application  to  alternative  missiles,  in 
addition  to  modifying  the  input  geometry  data,  the  trajectory  calculation  would 
need  to  be  modified.  There  are  a number  of  minor  changes  required  within  this 
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section  of  the  program;  for  instance  the  position  of  the  center  of  gravity  is 
defined  as  a linear  function  of  time  and  so  the  constants  involved  in  this 
expression  would  need  to  be  changed.  In  addition,  there  may  be  other  implicit 
assumptions  involved  within  the  code.  However,  the  main  change  necessary  would 
be  to  replace  the  dynamic  and  aerodynamic  data  tables  at  the  beginning  of  the 
Main  program  with  coefficients  appropriate  to  a new  missile.  It  should  also 
be  pointed  out  that,  since  the  program  still  relies  heavily  on  these  data 
tables  for  the  computation  of  the  steady  flow  effects, imprc  d trajectory 
calculations  could  be  achieved  by  replacing  these  data  tables  with  improved 
wind-tunnel  results  as  they  become  available. 

8.5  Overlay  Description 

To  reduce  the  computer  core  size  requirement  this  program  is  designed  to 
run  using  the  overlay  structure  defined  below.  When  run  in  this  way,  a region 
of  approximately  300K  bytes  of  computer  core  is  required. 
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8.6  Program  Test  Case 

A test  case  has  been  run  using  the  trajectory  input  data  described  in 
subsection  8.2  (with  HELIC=1.0)  together  with  the  following  panel  method  input 
data : 

(Card  Column) 

1 2 3 4 5 6 

1234 0123 0123 0123 0123 0123 012345 

0 0 

10  X 6 MISSILE  - ONE  PLANE  OF  SYMMETRY  TEST  CASE  2 

10  6 1.0  -1.0 

4.5  0.0  4.1625  0.0572915  3.8250  0.114583 

3.425  0.114583  2.9  0.114583  2.4  0.114583 

1.8  0.114583  1.2  0.114583  0.6  0.114583 

0.2  0.114583  0.0  0.114583 

This  represents  a 60  panel  test  case  with  a missile  length  of  4.5  feet  and 
diameter  of  2.75  inches.  For  production  runs,  the  precise  missile  geometry 
corresponding  to  the  2.75  Inch  Rocket  would  need  to  be  substituted. 

Some  of  the  output  is  given  in  figures  18(a)  and  (b).  Figure  18(a)  lists 
the  output  produced  by  the  panel  method  consisting  of  the  geometric  data 
defining  the  coordinates  of  the  panel  corners  together  with  the  components  of 
the  unit  normal  vectors  for  each  panel.  The  values  NPX,  NPY  and  NPZ  define 
the  control  point  coordinates  for  each  panel  while  the  final  column  lists 
some  other  geometric  properties  of  the  panels  including  their  area  which  are 
used  during  the  formation  of  the  velocity  matrices  within  the  panel  method. 

The  remainder  of  the  output  represents  a sample  of  the  output  produced 
during  the  trajectory  calculation,  of  which  the  first  page  and  the  last  two 
pages  have  been  included. 
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Figure  2.  Instantaneous  pitch  angle  for  triangular  downwash 


Instantaneous  pitch  angle  for  helicopter  downwash 
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Instantaneous  pitch  angle  for  helicopter  downwash  with  no  separation 
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helicopter  downwash  with  rotor  blades  initially 


Figure  10.  Instantaneous  pitch  angle  for  helicopter  downwash  with  rotor  blades  initially  at  90  . 
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downwash  with 


Instantaneous  pitch  angle  for  helicopter  downwash  using  0.0125  sec.  time  step 


Figure  16.  Instantaneous  pitch  angle  for  helicopter  downwash  using  0.005  sec.  time  step 
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Figure  18a.  Panel  method  output  for  test  case. 
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Figure  18a.  (continued) 
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Figure  18a.  (continued) 
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Figure  18a.  (concluded) 
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Figure  18b.  Trajectory  output  data  for  test  case. 
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Figure  18b.  (continued) 
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Figure  18b.  (concluded) 


